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Abstract 

Methylamine  radicals  (CH3NH)  and  amino  radicals  (NH2)  are  the  major  products  in  the  early 
ignition/pyrolysis  of  monomethylhydrazine  (CH3NHNH2).  The  kinetics  of  thermal  decomposition  of 
CH3NH  radicals  was  analyzed  by  RRKM  master  equation  analysis.  It  was  found  that  P  scission  of  the 
methyl  H  atom  from  CH3NH  radicals  is  dominant  and  fast  enough  to  induce  subsequent  H-abstraction 
reactions  to  trigger  ignition.  The  gas-phase  kinetics  of  H-abstraction  reactions  from  CH3NHNH2  by  H 
atoms  was  further  investigated  by  ab  initio  second-order  multireference  perturbation  theory  and 
quadratic-configuration-interaction  calculations.  It  was  found  that  the  energy  barriers  for  abstraction  of 
the  central  amine  H  atom,  two  terminal  amine  H  atoms,  and  methyl  H  atoms  are:  2.95,  4.16,  5.98,  and 
8.50  kcal  mof^,  respectively.  In  units  of  cm^  molecule'^  s'\  the  corresponding  rate  coefficients  were 
found  to  be:  kv  =  9.63xI0-^V-^‘^'’exp(-I54.2/T);  ks  =  2.04xI0''^T^'^^^  exp(I04.I/T);  kg  =  1. 13  x  I0'^° 
exp(-4I6.3/T);  kw  =  2.41x10'^^  j3.650  exp(-870.5/T),  respectively.  The  results  reveal  that 
abstraction  of  the  terminal  amine  H  atom  to  form  trans-CH3NHN»H  radicals  is  the  dominant  channel 
among  the  different  abstraction  channels.  At  298K,  the  total  theoretical  H-abstraction  rate  coefficient, 
calculated  with  no  adjustable  parameters,  is  8.16x10'  cm  molecule'  s'  ,  which  is  in  excellent 
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agreement  with  the  experimental  observation  of  (7.60  ±  1.14)xl0'  cm  molecule'  s'  (Vaghjiani,  J. 
Phys.  Chem.  A,  1997). 

Key  words:  CH3NH  decomposition,  CH3NHNH2,  H-abstraction,  CASPT2,  rate  coefficient 
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1,  Introduction 


The  binary  system  of  monomethylhydrazine  (CH3NHNH2)  and  dinitrogen  tertroxide  (N2O4)  is  widely 
used  as  liquid  propellants  for  satellite  and  spaee  propulsion  applieations  [1].  It  is  known  that  the  N — N 
bond  strength  is  weaker  than  that  of  the  C — N  bond  in  CH3NHNH2,  and  as  sueh,  the  CH3NH  and  NH2 
radieals  from  N — N  bond  fission  are  the  major  produets  in  the  early  ignition/pyrolysis  of  CH3NHNH2 
[2],  Using  the  high-pressure  limit  rate  eoeffieients  of  CH3NHNH2  decomposition  [3],  the  branching 
ratio  between  the  N — N  bond  fission  and  C — N  bond  fission  was  found  to  be  95:5  at  a  temperature  of 
lOOOK  and  for  pressures  in  the  range  of  2  to  3  atm.  To  experimentally  quantify  the  pressure-dependent 
rate  coefficients  for  the  N — N  bond  fission,  Hanson’s  group  [4]  recently  measured  the  NH2 
concentration  time-histories  in  CH3NHNH2/Ar  mixtures  over  the  temperature  range  of  1100  to  1400K 
and  pressures  in  the  range  of  0.3  to  5  atm  behind  reflected  shock  waves  using  laser  absorption 
measurement  methods.  The  fate  of  the  CH3NH  radicals  however  was  not  reported.  Theoretically,  we 
have  found  that  the  CH3NH  radical  is  very  reactive  since  direct  P-scission  of  the  methyl  H  atom  in 
CH3NH  has  a  small  energy  barrier  to  form  amino-methylene  (CH2=NH)  as  the  stable  product. 
Furthermore,  the  H  atoms  generated  from  this  decomposition  reaction  will  react  with  the  unbumed  fuel 
molecules  to  induce  ignition. 

In  this  work,  the  kinetics  of  the  thermal  decomposition  of  CH3NH  radicals  was  analyzed  and 
pressure  dependent  rate  coefficients  for  the  dominant  reaction  pathways  determined  by  RRKM  master 
equation  analysis.  Furthermore,  the  kinetics  of  H-abstraction  from  CH3NHNH2  by  H  atoms  was 
analyzed  using  second-order  multireference  perturbation  theory  and  unrestricted  quadratic- 
configuration- interaction  QCISD(T)  calculations.  The  rate  coefficients  for  each  of  the  abstraction 
channels  were  quantified  by  micro-canonical  transition  state  theory,  and  the  associated  branching  ratios 
determined.  It  was  found  that  the  total  rate  coefficient  determined  from  ab  initio  kinetics  theories  is  in 
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excellent  agreement  with  the  experimental  value  observed  at  room  temperature  [5].  Thus  the  theoretical 
rate  coefficients  from  this  work  should  be  reliable  over  a  wide  temperature  range  and  are  essential 
parameters  in  the  development  of  reaction  mechanisms  for  the  kinetics  modeling  of  ignition  involving 
monomethylhydrazine  and  its  derivatives. 

2,  Electronic  Structures 

There  are  several  reaction  pathways  for  the  CH3NH  radical  decomposition,  which  include  P 
scission,  isomerization,  molecular  elimination,  and  bond  fission.  The  geometries  of  the  stationary  points 
of  the  potential  energy  surfaces  (PESs)  for  these  reactions  were  optimized  at  the  B3LYP/6-31 1++G(d,p) 
level,  except  for  bond  fission  reactions.  For  the  C — N  and  N — H  bond  fission  reactions  in  CH3NH, 
multi-reference  character  of  the  wave  functions  is  significant.  Hence  geometries  of  the  transition  states 
for  these  two  bond  breaking  channels  were  optimized  by  multireference  second-order  perturbation 
theory  [6]  with  appropriate  orbitals  included  in  the  active  space.  The  C — N  and  N — H  bond  fissions 
involve  doublet  and  quartet  PESs  because  of  the  parallel  and  antiparallel  alignment  of  electron  spins  in 
the  NH  and  CH3N  radicals.  Since  the  energy  of  NH(X  )  is  42.73  kcal/mol  lower  than  that  of 
NH(a  A),  the  reaction  path  leading  to  the  CH3  +  NH(X  S  )  products  is  energetically  favorable.  This  is 
also  true  in  the  case  for  the  triplet,  CH3N  product  formation.  Therefore,  the  transition  state  structures  in 
C — N  and  N — H  bond  fission  reactions  were  determined  for  the  triplet  products  formation.  Specifically, 

-5  _ 

for  the  CH3  +  NH(X  E  )  dissociation  channel,  the  active  space  includes  the  bonding  orbital  between  the 
p  orbital  of  the  N  atom  and  the  p  orbital  of  the  C  atom,  as  well  as  a  perpendicular  p  orbital  for  the  radical 
on  the  NH  fragment.  Similarly,  for  the  CH3N  +  H  dissociation  channel,  the  active  space  includes  the 
bonding  orbital  between  the  p  orbital  of  the  N  atom  and  the  s  orbital  of  the  H  atom,  as  well  as  a 
perpendicular  p  orbital  for  the  radical  on  the  CH3N  fragment.  Such  orbital  configurations  were  included 
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in  the  active  space  of  the  transition  state  to  acquire  a  parallel  alignment  of  electron  spins  on  the 
corresponding  N  atom  for  the  formation  of  the  triplet  NH  and  CH3N  products. 

For  the  abstraction  of  amine  hydrogen  atoms  by  H  atom,  the  multi-reference  character  of  the 
wave  functions  was  also  found  to  be  quite  significant,  as  reflected  by  their  T1  diagnostics  [6]  of  0.03 
from  QCISD(T)/cc-pVQZ  calculations.  The  CASPT2  method  with  aug-cc-pVDZ  and  aug-cc-pVTZ 
basis  sets  was  applied  to  optimize  the  geometries  of  reactants  and  transition  states,  and  to  calculate  the 
ro-vibrational  frequencies.  In  the  CASPT2  geometry  optimizations,  the  active  space  (4e,3o)  was  chosen 
for  CH3NHNH2,  consisting  of  the  a,  g*  orbital  pair  of  the  N — N  bond  and  the  s  orbital  of  the  amine  N 
atom.  The  active  space  (3e,3o)  consisting  of  the  H  atom  orbital  and  the  a,  a*  orbital  pair  of  the  H  atom 
being  abstracted  was  chosen  for  the  transition  states  for  the  abstraction  of  amine  H  atoms.  For 
abstraction  of  the  methyl  H  atom,  the  orbitals  of  the  N — N  and  C — N  bonds  become  delocalized. 
Consequently,  a  larger  active  space  (7e,5o)  consisting  of  the  H  atom  orbital,  the  a,  a*  orbital  pair  of  the 
methyl  H  atom  being  abstracted,  and  the  two  hybrid  orbitals  between  the  methyl  and  the  NH  groups  was 
chosen  for  the  transition  state  for  the  abstraction  of  methyl  hydrogen. 

Higher-level  single-point  energies  were  determined  from  the  spin  unrestricted  or  restricted 
quadratic-configuration-interaction  QCISD(T)  calculations  for  open-shell  and  closed-shell  species.  The 
cc-pVTZ  and  cc-pVQZ  basis  sets  were  adopted  in  the  QCISD(T)  calculations  and  the  energies  were 
extrapolated  to  the  complete  basis  set  (CBS)  limit,  cc-pVooZ,  by  the  asymptotic  form  [7]  [8].  Electronic 
structure  calculations  were  performed  using  the  MOLPRO  [9]  and  Gaussian09  [10]  program  packages. 

3,  Kinetics  of  CH3NH  Radical  Decomposition 

The  thermal  decomposition  of  the  CH3NH  radicals  was  found  to  occur  by  the  following  reaction 
pathways: 

CH3NH  ^  TSl  ^  CH3  +  NH  (Rl) 
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^  TS2  ^  CH2=NH  +  H 

(R2) 

^  TS3  ^  CH2NH2 

(R3) 

CH2NH2 

^  TS4  ^  CH2=NH  +  H 

(R4) 

CH3NH 

^TS5  ^CH3N+H 

(R5) 

^  TS6  ^  CH2=N+  H2 

(R6) 

Figure  1  shows  the  stationary  points  of  the  potential  energy  surfaces  of  the  CH3NH  radical 
decomposition  characterized  at  the  UQCISD(T)//cc-pVooZ  level  with  CASPT2/cc-pVTZ  and  B3LYP/6- 
31 1++G(d,p)  geometries.  For  the  C — N  bond  fission  (Rl),  a  saddle  point  with  a  doublet  state  (TSl)  was 
found  to  have  an  energy  of  79.51  kcal/mol  higher  than  that  of  the  CH3NH  radical,  and  it  is  2.98  kcaFmol 
higher  than  that  of  the  CH3  +  NH(X^S~)  channel.  The  P  scission  of  the  methyl  H  atom  from  the  CH3NH 
radical  was  found  to  have  an  energy  barrier  of  33.08  kcal/mol,  forming  amino-methylene  (CH2=NH) 
(R2).  Isomerization  of  the  CH3NH  radical  produces  the  CH2NH2  radical  with  an  energy  barrier  of  36.39 
kcal  moF^  (R3).  The  CH2NH2  radical  formed  also  undergoes  P  scission  of  the  amine  H  atom  with  an 
energy  barrier  of  39.37  kcal/mol,  forming  CH2=NH  (R4).  The  N — H  bond  fission  reaction  in  the 
CH3NH  radical  that  leads  to  the  triplet  CH3N  radical  and  H  atom  is  endothermic  by  83.86  kcal  moF 
(R5).  The  doublet  transition  state  (TS5)  was  found  to  have  an  energy  of  0.54  kcal/mol  lower  than  its 
products.  A  direct  H2  elimination  from  the  CH3NH  radical  is  endothermic  by  13.24  kcal  mof'  (R6),  but 
its  energy  barrier  is  quite  high,  99.62  kcal  moF\  and  is  not  kinetically  relevant. 

In  the  decomposition  of  CH3NH  radicals,  micro-canonical  rate  coefficients  as  functions  of  total 
rotational-vibrational  energy  {E)  and  total  angular  momentum  (J)  were  calculated  using  Rice- 
Ramsperger-Kassel-Marcus  (RRKM)  theory  and  multi-well  master  equation  analysis  [11][12][13] 
which  were  implemented  in  the  VARIFLEX  code  [14].  Tunneling  corrections  were  included  for  all 

transition  states  based  on  asymmetric  Eckart  potentials  [15].  Pressure-dependent  kinetics  analysis  was 
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performed  by  solving  the  time-dependent,  multiple-well  master  equation.  The  collisional  energy  transfer 
probability  in  the  master  equation  analysis  was  approximated  by  a  single-exponential-down  model, 
employing  the  temperature  dependent  form,  AEdown  =  200(1/300)°'^^  cm'\  for  the  average  downward 
energy  transfer  [16].  The  Lennard- Jones  parameters  for  the  collision  rates  were  estimated  to  be  a  =  3.8 
A  and  s  =  253  cm'^  for  the  CH3NH  radical,  and  a  =  3.70  A  and  s  =  66.1  cm"'  for  nitrogen  [17]. 

We  chose  CH3  +  NH(X  CH3NH  (R-1)  as  an  entrance  channel  to  determine  the 

phenomenological  rate  coefficients.  At  a  pressure  of  1  atm,  the  rate  coefficients  for  three  dominant 
decomposition  channels  (R2,  R3,  and  R4)  were  obtained  and  fitted  as  the  following  expressions  in  the 
temperature  range  of  600  to  2500K. 

k2  =  2.69x10^’  exp(-18598.9/T)  s'^ 

k3=  1.02x10^°  t“  exp(-18823.9/T)  s'^ 

^4=  1.85x10^^  exp(-22570.3/T)  s'^ 

It  is  noted  that  due  to  the  uncertainties  in  the  parameters  for  the  average  energy  transfer  model,  the 
above  pressure  dependent  rate  coefficients  need  to  be  further  quantified/validated  by  future  experimental 
measurements.  Nevertheless,  these  rate  coefficients  should  be  useful  for  kinetics  analysis  for  the 
CH3NH  radical  decomposition.  Figure  2  shows  the  rate  coefficients  k2,  ks,  and  k4  at  1  atm  along  with  the 
corresponding  high-pressure  limit  rate  coefficients  in  the  range  of  600  -  2500K.  It  is  seen  that  the  direct 
dissociation  rate  (^2)  to  CH2=NH  +  H  is  dominant  at  all  temperatures.  The  dissociation  rate  (kA  to 
CH2=NH  +  H,  is  comparable  to  k2  only  when  the  temperature  is  above  '-'1800K.  The  isomerization  rate 
coefficient  {ks)  to  CH2NH2  is  about  one  order  of  magnitude  smaller  than  k2.  This  is  ascribed  to  the  tight 
three-center  transition  state  (TS3)  with  an  imaginary  frequency  of  1991.6  cm'\  while  the  P  scission  of 
the  methyl  H  atom  occurs  via  a  looser  transition  state  (TS2)  with  an  imaginary  frequency  of  555.4  cm'\ 
Furthermore,  strong  pressure  dependence  in  the  dissociation  rates  is  observed,  and  scission  of  the 
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methyl  H  atom  to  form  CH2=NH  dominates  at  all  pressures.  Between  1  atm  and  the  high  pressure  limit, 
the  rate  coefficient  k2  was  found  to  be  in  the  range  from  1.6x10  to  8.7x10  s'  ,  and  is  fast  enough  to 
produce  H  atoms  that  can  subsequently  abstract  H  atoms  from  unbumed  fuel  molecules  to  trigger 
ignition.  In  the  following  section,  the  kinetics  of  H-abstraction  from  CH3NHNH2  by  H  atoms  is 
analyzed. 

4,  Kinetics  of  H- Abstraction  from  CH3NHNH2  by  H  atom 

For  the  direct  abstraction  of  different  H  atoms  from  CH3NHNH2  by  H  atoms,  four  distinct  radicals  are 
possible,  as  shown  in  Figure  3. 

CH3NHNH2  +  H  ^  TS7  ^  CH3N*NH2  +  H2  (R7) 

^  TS8  ^  trans-CH3NHN*H  +  H2  (R8) 

^  TS9  ^  cis-CH3NHN*H  +  H2  (R9) 

^TS10^C*H2NHNH2  +  H2  (RIO) 

For  the  CH3NHN*H  radical,  it  is  noted  that  the  geometry  with  the  terminal  amine  H  atom  trans 
to  the  methyl  group  is  defined  as  the  trans-CH3NHN*H  radical,  and  the  geometry  with  the  terminal 
amine  H  atom  cis  to  the  methyl  group  is  defined  as  the  cis-CH3NHN*H  radical.  In  these  CASPT2 
optimized  transition  state  geometries,  the  cleaving  N — bond  lengths  are  1.137,  1.128,  and  1.163  A, 
and  the  forming  H — bond  lengths  are  1.071,  1.127,  and  1.038  A  in  TS7,  TS8,  and  TS9,  respectively. 
In  TSIO,  the  cleaving  C — bond  length  is  1.232  A  and  the  forming  H — bond  length  is  1.172  A.  In 
comparing  the  above  cleaving  and  forming  bond  length  data,  it  is  seen  that  the  formation  of  the  trans- 
CH3NHN*H  radical  via  TS8  should  be  more  facile  because  of  the  longer  forming  H — bond  length. 
This  implies  a  low  energy  barrier  in  this  channel  and  as  such  is  the  dominant  reaction  path. 
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With  the  optimized  geometries  in  Figure  3,  the  rotational  and  vibrational  frequeneies  of  the  four 
transition  states  were  ealculated  at  the  same  CASPT2/ec-pVTZ  level.  The  T1  diagnosties,  zero-point 
energies,  and  zero-point  eorreeted  UQCISD(T)/co-pVooZ  energies  relative  to  that  of  CH3NHNH2  +  H 
ehannel  are  listed  in  Table  1.  The  four  abstraetion  channels  (R7  through  RIO)  have  energy  barriers  of 
4.16,  2.95,  5.98,  and  8.50  kcal  mof',  respectively.  The  channels  for  the  abstraction  of  amine  H  atoms 
are  exothermic  by  23  to  25  kcal  mof^  while  the  channel  for  the  abstraction  of  the  methyl  H  atom  is  less 
exothermic  (11.7  kcal  mol'^).  The  abstraction  of  a  terminal  amine  H  atom  to  form  trans-CHsNHN'H  + 
H2  (R8)  has  the  lowest  energy  barrier  of  2.95  kcal  mof'  among  the  four  transition  states,  which  is  due  to 
the  longer  forming  H — bond  in  TS8  as  discussed  above.  It  is  expected  to  be  the  dominant  channel 
among  these  abstraction  reactions. 

The  rate  coefficients  of  the  above  H-abstraction  reactions  were  determined  using  {E,  J)  resolved 
micro-canonical  transition  state  theory  with  the  CASPT2/aug-cc-pVTZ  optimized  geometries  and 
vibrational  frequencies.  Hindered  rotor  corrections  for  the  C-N  and  N-N  torsional  modes  were 
obtained  from  one-dimensional  fits  of  the  torsional  potentials,  employing  Pitzer-Gwinn  like 
approximations  [18]  and  the  moments  of  inertia  that  reproduced  the  coupled-harmonic-oscillator 
limit  at  low  temperatures,  and  the  free-rotor  limit  at  high  temperatures.  Tunneling  effects  through  the 
saddle  points  were  modeled  with  asymmetric  Eckart  potentials.  The  imaginary  frequencies  of  the  four 
transition  states  used  for  the  tunneling  corrections  are:  1895.4,  1674.7,  1969.9,  and  2151.2  cm'\ 
respectively,  which  were  obtained  from  the  CASPT2/cc-pVTZ  calculations.  The  calculated  rate 
coefficients  for  the  abstraction  of  different  H  atoms  in  CH3NHNH2  by  H  atom  are  shown  in  Figure  4. 
The  rate  coefficient  kg  for  the  abstraction  of  the  terminal  amine  H  atom  to  form  trans-CH3NHN*H  +  H2 
dominates  below  1500K.  Furthermore,  the  rate  coefficient  kg  is  within  the  lower  limit  of  the  reported 
experimental  rate  coefficient  at  298K  [5],  which  implies  that  the  channel  (R8)  is  the  predominant 
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abstraction  channel.  Due  to  the  small  energy  barrier,  tunneling  effeets  on  k%  are  expected  to  be 
significant  at  low  temperatures.  The  tunneling  faetor  for  was  found  to  be  in  the  range  of  6.20  -  1.71 
for  temperatures  in  the  range  of  298  -  600K,  and  deereases  to  1.37  -  1.04  in  the  higher  temperature  range 
of  800  -  2500K.  The  abstraction  of  the  central  amine  H  atom  (R7)  beeomes  dominant  above  1500K, 
whereas  reaetions  (R9)  and  (RIO)  become  important  at  temperatures  above  lOOOK.  It  is  seen  that  the 
reactions  are  governed  by  the  abstraction  of  the  amine  hydrogen  atoms  at  low  temperatures,  and  the 
abstraetion  of  the  methyl  H  atoms  becomes  important  at  high  temperatures.  Furthermore,  P-scission  of 
NH2  from  the  C*H2NHNH2  radieal  is  the  fastest  proeess  among  the  decomposition  channels  of  the  four 
radieals  [19].  Therefore  all  of  the  abstraetion  channels  are  important  pathways  in  the  CH3NHNH2 
pyrolysis/oxidation.  In  units  of  cm  molecule'  s'  ,  the  rate  coefficients  were  fitted  as  follows: 

^7=  9.63  X  10-20  exp(-154.2/T) 

k8  =  2M  X  exp(104.1/T) 

k9=  1.13  X  exp(-416.3/T) 

kio  =  2.41  X  10'^^  exp(-870.5/T) 

Assuming  the  rate  coefficients  are  additive,  the  total  rate  coefficient  ktotai  was  determined  to  be 
2.55  X  10'^°  T^'^^^exp(227.6/T)  em^  moleeule'^  s'\  At  room  temperature,  the  theoretical  rate  eoefficient 
IS  8.16  X  10'  cm  molecule'  s'  ,  and  as  seen  in  Figure  4,  it  is  in  the  upper  bound  of  the  experimental 
observation  of  (7.60  ±  1.14)  X 10'  cm  molecule'  s'  [5].  The  excellent  agreement  between  theory  and 
experiment  for  the  total  rate  coefficient  at  room  temperature  gives  us  eonfidence  in  the  aceuracy  of  the 
theoretical  rate  eoefficients  at  high  temperatures.  Based  on  these  theoretieal  rate  coefficients,  the 
branehing  ratios  as  a  funetion  of  temperatures  are  plotted  in  Figure  5.  k^,  for  forming  the  trans- 
CH3NHN*H  radieal  is  shown  to  be  the  predominant  rate  coefficient  among  all  of  the  abstraction  channel 


DISTRIBUTION  A:  Approved  for  public  release,  distribution  unlimited 


rate  coefficients  for  temperatures  below  1200K.  This  is  consistent  with  our  previous  study  on  the 
related  H  atom  abstraction  from  CH3NHNH2  by  OH  radicals  [20]. 

Recently,  a  kinetics  analysis  on  the  reactions  of  CH3NHNH2  +  H/D  was  performed  by  Wang  et 
al.  [21]  based  on  quantum  chemistry  calculations  at  the  MCG3-MPWPW91//MPW1K/6-31  lG(d,p) 
level,  which  yielded  a  value  for  the  total  rate  coefficient  of  1.47xlO'^^T^  ^"^exp(-239.1/T)  cm^  molecule"' 
s"'  in  the  temperature  range  of  200  to  2000K.  Their  calculated  total  rate  coefficient  at  298K  is  4.14  xlO" 
cm^  molecule"'  s"',  which  is  a  factor  of  2  lower  than  the  experimental  value  [5].  The  discrepancy  in 
the  total  rate  coefficient  value  between  the  study  of  Wang  et  al.  [21]  and  this  work  is  in  the  range  2  -  0  .7 
in  the  temperature  range  of  300  -  2500K.  In  their  analysis,  Wang  et  al.  considered  the  H-abstraction 
reactions  of  two  CH3NHNH2  isomers  with  six  distinguishable  pathways  for  each  isomer.  The  rate 
coefficient  for  each  reaction  channel  was  calculated  by  using  canonical  variational  transition  state  theory 
with  a  small-curvature  tunneling  correction  (CVT/SCT),  and  all  of  the  rotational  and  vibrational 
frequencies  were  treated  by  the  harmonic-oscillator  approximation.  The  total  rate  coefficient  was 
determined  by  summation  of  the  individual  rate  coefficient  contributions  from  each  of  the  isomers  using 
weighting  factors  in  the  Boltzmann  distribution  function.  Furthermore,  they  concluded  that  the  H  atom 
abstraction  from  the  central  NH  group  is  the  dominant  channel  below  500K,  which  is  contrary  to  our 
findings.  We  note  here  that  the  branching  ratios  will  certainly  affect  the  radical  product  distributions 
when  modeling  CH3NHNH2  pyrolysis  and  oxidation. 

In  this  work,  we  have  carried  out  high-level  UQCISD(T)/cc-pVooZ//CASPT2/aug-cc-pVTZ 
calculations  to  obtain  accurate  values  for  the  molecular  properties  and  energy  values  for  the  species 
involved  in  this  system.  The  rate  coefficient  was  computed  by  micro-canonical  transition  state  theory 
with  conservation  of  total  rotational-vibrational  energy  {E)  and  total  angular  momentum  (J),  and  the 
torsional  anharmonicities  were  treated  as  internal  rotations  for  the  rotational  conformers.  Consequently, 
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our  approach  should  provide  a  better  perspective  for  the  branching  ratios  as  well  as  absolute  values  of 
the  temperature  dependent  rate  eoefficients. 


5,  Concluding  Remarks 

The  kinetics  of  the  thermal  deeomposition  of  CH3NH  radicals  was  analyzed  and  pressure  dependent  rate 
coeffieients  of  the  dominant  reaetion  pathways  were  determined  by  RRKM  master  equation  analysis.  P 
seission  of  methyl  H  atom  to  form  the  CH2=NH  product  was  found  to  be  the  most  dominant  channel  in 
the  CH3NH  radical  decomposition.  Furthermore,  the  kineties  of  H-abstraction  from  CH3NHNH2  by  H 
atoms  was  analyzed  using  seeond-order  multireferenee  perturbation  theory  and  unrestricted  quadratie- 
configuration-interaetion  QCISD(T)  ealculations.  The  rate  coefficient  for  eaeh  ehannel  was  calculated 
at  the  {E,  J)  resolved  level  using  micro-canonical  transition  state  theory,  and  the  corresponding 
branehing  ratios  for  the  produet  channels  determined.  The  results  reveal  that  the  abstraetion  of  the 
terminal  amine  H  atom  to  form  trans-CH3NHN*H  radical  is  the  most  dominant  channel  among  all  the 
abstraetion  ehannels.  The  total  theoretical  rate  coefficient  eomputed  is  in  exeellent  agreement  with  the 
available  experimental  data. 
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TABLE  1  ■  Molecular  Properties  for  CH3NHNH2  +  H  Reaction 


Transition  state 

T1  diagnostic  ^ 

zpe'’ 

AE" 

AE‘^ 

TS7_CH3N*NH2 

0.029 

49.76 

4.16 

4.27 

TS8  _trans-CH3NHN*H 

0.031 

50.13 

2.95 

3.44 

TS9_eis-CH3NHN*H 

0.028 

49.79 

5.98 

5.91 

TSIO  C-H2NHNH2 

0.018 

50.08 

8.50 

8.11 

T1  diagnostics  calculated  at  the  QCISD(T)/cc-pVQZ  level.  Zero-point  energy  ealeulated  at  the 
CASPT2/aug-ec-pVTZ  level,  units  in  kcal  mol’'.  Zero-point  eorreeted  energy  relative  to  the  reaetants 
CH3NHNH2  +  H  ealeulated  at  the  QCISD(T)/cc-pVooZ//CASPT2/aug-ee-pVDZ  level,  units  in  keal  mol' 
Zero-point  eorreeted  energy  relative  to  the  reaetants  CH3NHNH2  +  H  ealeulated  at  the  MCG3- 
MPWPW91  //MPWlK/6-31  lG(d,p)  level  by  Wang  et  al.  [21],  units  in  keal  moP'. 
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Figure  5 


Figure  Captions 

Figure  1.  Potential  energy  surfaees  for  CHaN'H  radieal  decomposition. 

Figure  2.  Rate  coefficients  of  CH3NH  radical  dissociation  in  the  temperature  range  600  - 
2500K. 

Figure  3.  Potential  energy  surfaces  for  CH3NHH2  +  H  reaction  calculated  at  the 
UQCISD(T)/cc-pVooZ//CASPT2/aug-cc-pVTZ  level. 


Figure  4.  Theoretical  abstraction  rate  coefficients  for  CH3NHH2  +  H  reaction, with  comparison 
of  experimental  data. 

Figure  5.  Theoretical  abstraction  branching  ratios  for  CH3NHH2  +  H  reaction  as  a  function  of 
temperature. 
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